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Abstract. We investigate critical wetting transitions for fluids adsorbed in wedge- 
like geometries where the substrate height varies as a power-law, z{x,y) ~ ja^T, in 
one direction. As 7 is increased from to 1 the substrate shape is smoothly changed 
from a planar-wall to a linear wedge. The continuous wetting and filling transitions 
pertinent to these limiting geometries are known to have distinct phase boundaries 
and critical singularities. We predict that the intermediate critical wetting behaviour 
occurring for < 7 < 1 falls into one of three possible regimes depending on the 
values of 7,p and q. The unbinding behaviour is characterised by a high degree of 
non-universality, strongly anisotropic correlations and enhanced interfacial roughness. 
The shift in phase boundary and emergence of universal critical behaviour in the linear 
wedge limit is discussed in detail. 
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There is growing interest from theorists and experimentahsts in fluid adsorption on 
micro-patterned and sculpted solid substrates |l|, ^, H, ||, ^ ||, ^ |TT|, |12|, 
|14|, |T^. Surface decoration and structure can substantially alter the character of fluid 
adsorption and lead to novel examples of interfacial phase transitions and enhanced 
fluctuation related effects. There is also strong evidence that there are fundamental 
connections between geometry-induced and fluctuation-induced interfacial phenomena 
at wedge filling transitions [0] and also between wedge filling and unzipping transitions 



for doubled-stranded DNA ig. 

The purpose of the present article is to focus on interfacial adsorption in generalised 
3D wedge shaped geometries, which are translationally invariant in one direction (along 
the y axis, say), but whose height above some reference plane varies as a power law 
z{x,y) ~ for large distances in the x direction. We refer to this particular 

class of surface geometry as the gamma-wall which may be viewed as an example of 

f To whom correspondence should be addressed (anna.sartori@curie.fr) 
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deterministic roughness |T7| , |T^, |T9|. Note that by increasing the exponent 7 the wall 
morphology can be changed smoothly from a planar substrate (7 = 0) to a linear 
wedge (7 = 1) and eventually to a parallel plate geometry (7 = 00). The adsorption 
properties and interfacial fluctuation effects in each of these geometries, corresponding 
to wetting |^ filling |1^, |1T| and capillary condensation respectively, are 
very different to each other and have received considerable theoretical and experimental 
interest. The central question we ask here is, how do the wetting properties depend on 
the wall exponent 7? For substrates that are completely wet by the fluid (corresponding 
to vanishing contact angle 6 = 0), recent work [T^, |l^ has shown that the adsorption 
isotherms for the gamma wall show a sensitive dependence on 7, which facilitate the 
cross-over from continuous complete wetting (7 = 0) to first-order capillary condensation 
(7 = 00) phenomena, through a sequence of novel interfacial behaviours, which emerge 
at intermediate values of 7. Here, we extend this study to the case where the (planar) 
substrate undergoes a continuous (critical) wetting transition at some temperature T„ 
and ask how the phase boundary and critical exponents characteristic of the interfacial 
unbinding depend on 7. As we shall show this case is considerably more complex 
than the complete wetting scenario, due to the presence of large scale interfacial 
fluctuation effects, which gives rise to strongly non-universal critical behaviour and 
distinct fluctuation regimes. Nevertheless, the critical properties characteristic of these 
regimes are precisely those that allow us to understand the change in phase boundary 
and fluctuation-related properties, that occur as the wall morphology is changed from 
a planar substrate (7 = 0) to a linear wedge (7 = 1). 

Consider the interface between a non-planar substrate modelled as an inert 
spectator phase, whose height is described by the a continuous function z{x, y) in the 
shape of a generalised wedge (see figure |I|), with the power law behaviour described 
above. The substrate is in contact with a bulk vapour phase at temperature T and 
pressure p, which we will suppose is tuned to bulk two-phase coexistence p = psat(T). 
The planar wall-fluid interface (corresponding to 7 = 0) is taken to have a continuous 
wetting transition at temperature at which the contact angle 0{T) vanishes and the 
mean interfacial height diverges. The mean-field critical singularities occurring at the 
wetting transition are found by minimising the binding potential 

m) = -|; + ^. (1) 

where t oc (T^ — T)/T^, b > and the exponents p, q > p depend on the ranges of 
the intermolecular forces. The critical exponents describing the divergence of the mean 
interfacial height ~ t'^", roughness C,± ~ t~'^^ and parallel correlation length ~ 
are given by the well known expressions 

A = = o(v^), z/ii = (2) 

q-p V » 2{q-p) 

and are not altered by thermal fluctuation 2T|. Note that the interfacial roughness 
C,± is negligible compared to the wetting film thickness. Also recall that the contact 
angle vanishes as ^ ~ i(2-as)/2^ with 2 — = q/ {q — p). 



Letter to the Editor 



3 



Elementary thermodynamic considerations, based on balancing bulk and surface 
tension contributions to the grand potential, show that the location of the surface phase 
boundary/wetting transition in the gamma- wall geometry depends qualitatively on the 
behaviour of r = A-y/A^^, corresponding to the ratio of the total to planar (projected) 
area of the substrate. Thus, for < 7 < 1 for which r = 1, the wetting phase boundary 
remains at T = T^, whilst for 7 > 1, for which r = 00, the wedge is completely filled 
by liquid (gas) at all temperatures such that 0{T) is less (greater) than tt/2. Thus, for 
7 > 1 the wetting transition at bulk coexistence is superseded by a first-order capillary 
condensation-like phenomena. This is closely related to unbending phase transition 
occurring on corrugated surfaces |^. Hereafter, we restrict our attention to the regime 
7 < 1 and focus on the fate of the planar wetting transition as 7 is increased from to 
1. 

The critical behaviour occurring at the limit 7 = 1, corresponding to the filling of a linear 
wedge, is known in some detail p|, 0, |ll|]. Writing the wall- function z{x, y) = tana \x\, 
with a the tilt angle, observe that the ratio r = sec a > 1. Accordingly, the phase 
boundary for the wedge wetting (filling) is shifted and again thermodynamic arguments 
dictate that it occurs at a lower filling temperature Tf, satisfying 0{Ti) = a rather than 
6(T) = For planar substrates that undergo critical wetting transitions, the linear 
wedge filling transition is also continuous, but is characterised by critical exponents 



distinct to those at critical wetting |in, 11]. Consider, for example, the divergence 



of the wedge mid-point interfacial height /„ ~ t'~^"", mid-point roughness ^± ~ t'~'^^ 
and correlation length C,y ~ t''^^, where the latter is measured along the wedge and 
t' oc (Tf — T)/Tf. Calculations based on effective interfacial models show that the 
criticality falls into two regimes [l^. For p < 4 the critical exponents belong to a filling 
mean-field (FMF) regime with critical exponents /5„ = 1/p, = 1/4, Uy = 1/2 + 1/p, 
in which the roughness is much smaller than the film thickness, C,±/lw << 1- For P > 4, 
on the other hand, there is a filling fluctuation (FFL) regime and the critical exponents 
take the universal values 

(3^ = 1/4, = 1/4, Uy = 3/4. (3) 

In this regime, ~ /„ and the interfacial fluctuations are controlled by an effective 
wedge wandering exponent, so that /w ~ ^± ~ (,y with ^ = 1/3 . The universal value of 
z/^ implies that fluctuation effects, at the filling of a linear wedge, are always important 
and contrast sharply with those at planar wetting transition. 

With these preliminaries aside, we can now precise the two central questions 
addressed in this paper. First, we introduce critical exponents for the mid-point height, 
roughness and correlation length as t — >■ by the identifications 

L ~ t-f'^^^ e± ~ t-''^^'\ ~ t-^'^^\ (4) 

where we have restricted ourselves to the (unknown) critical behaviour occurring in 
the range < 7 < 1. Then, in evaluating the full 7, p and q dependence of these 
critical exponents, we wish to understand two specific points; A: How does the critical 
behaviour reflect the discontinuous shift of the phase boundary from 6' = to ^ = a as 
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7 ^ 1^ ? B: How do the critical exponents change from their wetting to filhng values? 
In particular, are any combinations of critical exponents continuous in this limit and 
allow us to smoothly 'turn on" the fluctuation effects? 

The starting point for our calculations is the interfacial Hamiltonian model 



H[l] = J Jdxdy 



^(Wlf + Wil-zix,y)) 



(5) 



where l{x,y) is measured relative to the horizontal reference plane ^ = 0, S denotes 
the stiffness (surface tension) of the unbinding liquid-vapour interface and W{1) is the 
binding potential (|I]). The model is only valid for substrate height functions z{x, y) that 
have a shallow gradient |V-2| << 1 and thus suffices to determine the critical behaviour 
in the regime of interest, < 7 < 1. The critical exponents are insensitive to the precise 
nature of the substrate shape near the x = line and the power-law shape may be cut 
off at some appropriate short- distance. 

The results of our analytical and numerical studies show that the critical 
behaviour is strongly non-universal and depends sensitively on the values of 7,p and 
q. We concentrate on the interpretation of these results, making brief mention of our 
calculational details at the end of our article. Figure ^ shows how the critical 
behaviour falls into three possible categories labelled the planar mean- field (IIMF), 
geometrical mean-field (GMF) and geometrical fluctuation (GFL) regimes, respectively. 
Within the IIMF regime, < 7 < 71(g) with 

7,(,) - j^. (6) 

the substrate shape does not alter the values of the wetting critical exponents and /5(7), 
z/_l(7) and 2^3/(7) are identical to /3s, z/_|_ and i/y shown in (|]). For 7 > 71, corresponding to 
the geometrical, GMF and GFL regimes, on the other hand, the wedge geometry alters 
the critical exponents from their planar values. In both these regimes the interfacial 
height diverges with a modified exponent 

Note that the critical exponent j3{'~f) > (3s, for 7 > 71, so that in the limit t ^ 
the mid-point height ^ In- In figure ^ we show numerical results for the mid-point 
height for the case of non-retarded van der Waals forces (p = 2, g = 3) and 7 = 3/5, 
for which we predict /3(3/5) = 9/4. Notice that the initial divergence of the interfacial 
height is planar-like (with critical exponent /5s = 1), but crosses over to the asymptotic 
geometrical result as t — * 0. The boundary between the GMF and IIMF regimes occurs 
when /3(7i) = /3s, so there is smooth cross-over in the critical behaviour for the interfacial 
height at the HMF/GMF separatrix. The GMF and GFL regimes are distinguished from 
each other by the behaviour of the correlation length critical exponents i^±(7) and 2^3/(7). 
For p < 4 the GMF regime extends from 71 to the linear wedge limit 7 = 1~. For p > 4:, 
on the other hand, the GMF regime terminates at 7 = 72, where 
' , 2p 
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Within the GMF regime, the roughness is still small compared to the film thickness, 
but is larger than in the IIMF regime and z/^^^(7) > 0. There is pronounced non- 
universality in this regime with 



^r'ii) = \ + ^^^^ 



2(1-7) 
P 

7 



(9) 



4 4 

and 

^^^"(7) = I + Pil) ( 1 + ) . (10) 



There is a smooth change from the IIMF to the GMF regime for the fluctuation critical 
exponents so that, at the separatrix, z/j^^^(7i) = 0. For 7 > 72 (relevant for systems 
with p > A only), corresponding to the GFL regime, on the other hand, the wedge 
wetting transition in the gamma wall is fluctuation dominated and we can identify 

.f^^ = /5(7), z^r^ = C(7)-r> (11) 

where (^(7) = 7/(7 + 2) is the wedge wandering exponent for the gamma wall. Thus, 
in the GFL regime one has simple scaling relations between the diverging lengthscales, 
~ ^± ~ ^y^'^^ ■ For systems with purely short-ranged forces, the GFL regime spans 
the entire range < 7 < 1 and /?(7) = z/y 7/ (1 — 7). 

Returning to the more general case observe again there is smooth cross-over between 
the behaviour of the fluctuation critical exponents at the separatrix between the GMF 
and GFL regimes with, for example, i^x^^(72) = /^(72)- The geometrical regimes are 
also characterised by a strong degree of anisotropy, with the correlation length C,y much 
greater than the lateral extent of the filled region ^3. ~ l^"'. Interfacial fluctuations 
within the GMF and GFL regimes are pseudo-one dimensional, in contrast with IIMF 
regime. 

At this point, a number of remarks are in order. 

(I) The existence of three fluctuation regimes for critical wetting in the generalised 



wedge geometry contrasts with the case of complete wetting 0, |T3l , for which there are 
only two and no significant enhancement of fluctuation effects. Intriguingly, the number 
of regimes for critical and complete wetting in the gamma-wall wedge are the same as 
that induced by thermal (or impurity induced) fluctuation effects at planar critical and 
complete wetting transitions, respectively [plj] . 

(II) In the geometry affected regimes (GMF and GFL), the equilibrium profile 
/eq(a;) has a particularly simple structure which can be seen to directly lead to the 
critical exponent identification (0). Near the center of the wedge, the interface is flat 
and at near constant height /eq( 

x) ~ /w At some distance ~ ll/'^, the interface 
strikes the wall and thereafter closely follows its shape. The interfacial height at the 
wedge mid-point is determined by the simple condition that the local angle of incidence 
between the interface and wall is equal to the contact angle. Observe that, as 7 — *• 1~, 
the critical exponent for the interfacial height diverges. Including amplitude factors, 
we find that the mid-point height diverges as L ~ {U/ty^^\ with tw a non-universal 
constant. As 7 — > 1^ this implies that the height becomes macroscopic for all t < t^, 
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which represents a shift of the phase boundary from t = to t = t^. This is equivalent 
to the shift of the phase boundary from 6 = to 6 = a for hnear wedge fiUing and 
answers our first question A. 

(Ill) Some aspects of the interfacial fluctuations show a smooth change from 
wetting to filling-like behaviour, as 7 is increased, and allow us to give a quantitative 
answer to question B. This is most simply seen in the wedge wandering exponent 
C(7) = 7/(7 + 2) pertinent to the GFL regime, which generalises the linear wedge 
result ( = 1/3. Less obvious is the behaviour of the critical exponent ratio P{'j)/i^±{'j), 
which also recovers the linear wedge result, such that 

7^1- z^±(7) ^± 
where the RHS is equal to 4/p and 1 for p < 4 and p > 4, respectively. 

For systems with non-retarded van der Waals forces {i.e. with binding potential 
exponents p = 2, g = 3) we make the following predictions. The planar result pertinent 
to the standard critical wetting transition /3s = 1, z/_l = and z/|| = 5/2 are unchanged 
within a IIMF regime corresponding toO<7<2/5. For 1 > 7 > 2/5 the transition 
belongs to the GMF and the critical exponents are geometry sensitive. For example, at 
7 = 1/2 we predict 

mm = I Mm = \, = \. (13) 

The value of the roughness critical exponent v^{l/2) = 1/4 is significant, since it is 
independent of the value of q and is therefore also valid for tricritical wetting. This 
degree of universality is similar to the true universality of predicted for linear wedge 
filling. 

To finish, we make brief mention of the methods used in our calculations. In both 
the IIMF and GMF regimes the roughness is much smaller than the interfacial height 
and mean-field methods are appropriate. Numerical results obtained by minimising (|^) 
are complemented by analytical approaches following approximate solution to the Euler- 
Lagrange equation, based on standard variational methods. This is straightforward in 
the GMF regime, since the profile has a particularly simple structure. The results for the 
correlation length critical exponent Uy and roughness exponent z/^ were first obtained by 
solving the Ornstein-Zernike equation for the structure factor S{Q), corresponding to the 
Fourier transform of the mid-point height-height correlation function {l{yi)l{y2)) with 
respect to wave-vectors along the wedge. The critical behaviour in the GMF and GFL 
regimes can also be described using an effective one- dimensional model Hamiltonian 
which describes the energy-cost of constrained interfacial configurations in terms 
of the mid-point height l{y) = l{x = 0,y) only. The model can be derived from the 
underlying interfacial model (1) using standard methods, which have been previously 
applied to the linear wedge problem [|l^. The dimensional reduction, explicit in this 
method, is justified by the extreme anisotropy of fiuctuations in the GMF and GFL 
regimes as t — 0. 
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The reduced dimensional effective interfacial Hamiltonian has the form 

2 



dy 




(14) 



where cr is a non-universal constant (proportional to S) and V^{1) is the wedge binding 
potential. Minimization of Vy{l) identically recovers the mean-field expression for in 
the IIMF and GMF regimes. For large / this has the expansion 



where and d^ are non-universal constants. In the limit 7^1 , we find 



(15) 
and 



d^ ^ t/{p ~ 1), so that both V^{1) and (|T^) smoothly recover the linear wedge model 
considered in |1^. From the one- dimensional model it is straightforward to derive all 
the critical exponents quoted earlier using standard methods. 

In summary, we have investigated the geometry dependence of critical wetting 
exponents for fluids adsorbed in power-law wedges. Our results show that surface shape 
has both a stronger (and subtler) effect on critical wetting than complete wetting tran- 
sitions, with criticality falling into three possible regimes which facilitate the crossover 
from planar wetting to linear wedge filling. 
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Figure 2. Critical regimes for the gamma-wall wedge. The curves ci and C2 
represent the separatrixes 71(g) and 72(^7,^), obtained ior q = p + n, with n 
fixed. 
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Figure 3. Log-Log plot of the mid-point height vs i for 7 = 3/5 
and dispersion forces, showing cross-over from planar-hke (IIMF) to GMF 
behaviour, with asymptotic criticahty at /3(3/5) =9/4. 



